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The ground entanglement and thermal entanglement in quantum mixed spin chains consisting 
of two integer spins 1 and two half integer spins | arrayed as|-|-l-lina unit cell 
with antiferromagnetic nearest-neighbor couplings Ji(J2) between the spins of equal (different) 
magnitudes, are investigated by adopting the log-negativity. The ground entanglement transition 
found here is closely related with the valence bond state transition, and the thermal entanglement 
near the critical point is calculated and shown that two distinct behaviors exist in the nearest 
neighbor same kind of spins and different kind of spins, respectively. The potential application of 
our results on the quantum information processing is also discussed. 
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I. INTRODUCTION 



Entanglement, one of the most striking features of 
quantum mechanics Q , has been recognized as an impor- 
tant resource for quantum information processing 0, • 
Recently, the entanglement in the systems of condensed 
matter physics such as the Heisenberg models has cause 
much^attentionjfl EEM j IMSEM! iii 
EE 113, EE El, Hfl HH HI, HI, [H, HI, H^. Natural en- 
tanglement in the thermal equilibrium of the spin chains 
coupled by the exchange interaction has been found and 
it is shown that the temperature or the external magnetic 
field can enhance the pairwise entanglement of spins 
The relations between the ground state entanglement and 
quantum phase transition have also be revealed for the 
classes of spin- 1/2 chains [H, 0, ^J. All results con- 
cerning the entanglement in the Heisenberg spin chains 
showed that the entanglement can make a bridge across 
the condensed matter physics and quantum information 
theory. 

One of the most famous results in quantum spin chains 
is the Haldane conjecture which implies that the half- 
odd-integer antiferromagnetic Heisenberg chains should 
have gapless spectrum with algebraic decay of correla- 
tions at zero temperature and integer spin ones should be 
gapped with exponential decay of correlations. Recently, 
some authors have addressed whether the entanglement 
length is finite or infinite in the Heisenberg spin-1 chain 
which is a gapped quantum Hamiltonian and proven the 
existence of an infinite entanglement length as opposed 
to their finite correlation length 23]. Besides the s pin- 
1 chains, recent experimental achievements [2^, l29l l3Clj 
have stimulated the interest to study both analytically 
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and numerically the mixed quantum spin chains, such 
as the alternating quantum spin chain with S 1 = 1/2 
and S* 2 = 1 with antiferromagnetic nearest-neighbor ex- 
change interaction or more complicated hypothetical ID 
periodic structures of different spins with ferrimagnetic 
properties |3l| . The topology of spin arrangements in the 
ferrimagnetic chain plays an essential role on the varia- 
tion of energy gap associated with Haldane conjecture 
and on the magnetism of ground states. The ground 
state of ferrimagnetic chain with two kinds of different 
spins is either a spin singlet or ferrimagnetic due to the 
setup of the bipartite lattices following Marshall theorem 
[32|. By applying the Lieb-Schultz-Mattis theorem [33| . 
one can find whether the ferrimagnetic chains have a gap- 
less excitation though the theorem failed to predict the 
existence of a gapped excitation. In this Letter, we con- 
sider the ground state entanglement and thermal state 
entanglement in an interesting model which is called as 
quantum mixed spin 1/2-1/2-1-1 chain by making use of 
the quantum Monte Carlo method with a improved loop 
cluster algorithm |34j . Over the past years, the quantum 
Monte Carlo method with a improved loop cluster algo- 
rithm has proven it to be a efficient numerical tool to deal 
with the complicated coupled spin chains So it is 

very credible to numerically investigate the nonclassical 
properties such as the entanglement of the mixed spin 
chain by adopting this method. 

For first time, we investigate the ground state entangle- 
ment of the mixed spin chain up to 128 sites by exploring 
the log-negativity as an entanglement measure. Two dis- 
tinct behaviors of pairwise entanglement are found. We 
provide a evidence that the ground state entanglement 
between two 1/2 spins or between 1/2 spin and 1 spin 
are closely related to the so-called VBS phase transition. 
At the point of VBS phase transition, the ground state 
between two 1/2 spins can be transferred to two differ- 
ent kinds of spins i.e. the spin 1/2 and spin 1 or vice 
versa. Furthermore, the thermal state pairwise entan- 
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FIG. 1: Graphical representation of the Hamiltonian in 
Eq.(l), the black and white circles represent spins | and 1, 
respectively. 

glement is also studied and whether the thermal entan- 
glement can be improved by increasing the temperature 
or not is also addressed, which strongly depends on the 
parameter a = 3-ij J\. 

This paper is organized as follows: In Sec. II, we briefly 
outline the basic content of the quantum mixed spin 
chain and investigate the ground state entanglement 
between two nearest-neighbor spins | or between two 
nearest- neighbor spins h and 1. In Sec. Ill, we investigate 
the thermal state entanglement between two nearest- 
neighbor spins i or between two nearest-neighbor spins 
^ and 1. It is shown that the entanglement can be en- 
hanced by improve the temperature. In Sec. IV, there are 
some conclusive remarks. 



II. The ground state entanglement and the VBS 
phase transition of quantum mixed spin chain 

Now we consider the quantum antiferromagnetic mixed 
spin S 1 — S 1 — S 2 — S 2 chain displayed by Fig.l, whose 
Hamiltonian can be given by [35l | 

N/4 

H = y^(Jl5 f j„_ 3 ■ Sl n _ 2 + JlS\ n _i ' ■Sfn— 1 
n=l 

+ J\Sl n _i ■ S"l n + J2S\ n ■ sl n+1 ). (1) 

Here, the periodic boundary condition is adopted. In the 
case with S 1 = \ and S* 2 = 1, the model has been studied 
by quantum Monte Carlo ( QMC ) simulation [3j| and the 
nonlinear a model (NLSM) [371 l88| . The ground state of 
this model is nonmagnetic with gapped excitations. The 
gap varies as a function of the parameter a = J2/J1, and 
vanishes at a critical point a c , which implies a quantum 
phase transition between two different VBS state. Since 
the VBS states of the mixed spin chains actually contain 
an ensemble of spin singlet state which can be regards as 
the abundance quantum information resource. It is very 
desirable to study how the parameter a affects the en- 
tanglement of this model or what is the relation between 
the gapped excitation and the thermal entanglement. 

We study both the ground state entanglement and 
the thermal state entanglement of this system. The 
state of a quantum system at thermal equilibrium can 
be described by the density operator p = iexp(— /3H), 
where j3 = with ks the Boltzmann's constant. 



Z = Tre~^ ff is the partition function. The entanglement 
in the ground state or thermal state is called ground state 
entanglement or thermal entanglement, respectively. For 
antiferromagnetic mixed spin chains, the ground state is 
not degenerate. So the ground state can be directly ex- 
pressed by the density operator p = exp(— (3H) with 
[3 — > 00. By utilizing the SU(2) invariance of this system, 
the reduced density matrix /Z 1 ' 1 ) for two spins | or p( 1,2 ' 
for one spin \ and one spin 1 can be obtained as follows 
0: 

P (M) =g (M) |0,0)(0,0|+ IMXM, (2) 

'' i- i 

and 

t=-l/2 i=-3/2 

where | J, J z ) denotes the state of total spin J and z com- 
ponent J z , and i? 1 - 1 ' 1 ) and g*- 1,2 - 1 are the function of tem- 
perature. For characterizing the ground state pairwise 
entanglement and the thermal state pairwise entangle- 
ment, we adopt the log-negativity as the entanglement 
measure, which has been proven to be an operational 
good entanglement measure. The log- negativity N(p r ) 
for the two particle reduced density operator p r is de- 
fined as [3{| 

V( Pr ) = log 2 ||^|| ; (4) 

where p^ is the partial transpose of p r and || denotes 
the trace norm of p^, , which is the sum of the singular val- 
ues of p^ For the isotropic mixed spin chains with SU(2) 
symmetry, the log-negativity between any two spins can 
be proven to be directly related with the two-particle 
correlation function. Therefore, we can directly calcu- 
late the log-negativity by slightly changing the program 
for the QMC. In Fig. 2, we numerically calculate the log- 
negativity N^' 1 ' of the reduced density operator of two 
nearest-neighbor spin i in a unit cell of the ground state 
and the log-negativity N 1 - 1 ' 2 ^ of two nearest-neighbor dif- 
ferent spin i and 1. It is shown that jV' 1 ' 1 ) decreases from 
1 to zero as the parameter a increases from to 1. At 
the point of VBS phase transition labelled by a c ~ 0.768, 
the entanglement experience a sharply descent. While 
the entanglement between two nearest-neighbor different 
spins i and 1 characterized by log-negativity TV^ 1 ' 2 ) is 
zero for small values of a, and near the critical point it 
increases with a and experience a sharply ascent at the 
critical point. Unfortunately, we find that the ground 
state entanglement only exists between two spins in the 
same unit cell. In the numerical calculation of the log- 
negativity of the ground state, we have set (3 = 200 and 
V = 128. We conjecture that the descent or ascent of log- 
negativity near the critical point can be more sharp if the 
value of (3 is chosen very very large, though this can con- 
sume more and more computation time. From Fig. 2, one 



can imagine that the quantum mixed spin chain in zero 
temperature can be utilized as a very good entanglement 
switch by controlling the parameter a. Recently, Pachos 
and Plenio have shown that some kinds of atomic lattice 
can be use to simulation the Hamiltonian of various con- 
trollable Heisenberg spin chain |2jJ. It can motivate us 
to realize the quantum mixed spin chain in the system of 
atomic lattice. Then we easily obtain the quantum mixed 
spin chain with controllable parameter a. The details of 
suggestion is discussed elsewhere. 



III. The thermal state entanglement of quantum 
mixed spin chain 

In what follows, we turn to consider the thermal entan- 
glement of the mixed spin chain. In Fig. 3, we calculate 
the log-negativity A 7 ^ 1,1 ) characterizing the entanglement 
of two nearest-neighbor as the function of ksT. It is 
shown that, if the parameter a is smaller than the crit- 
ical point a c , the entanglement decreases with the tem- 
perature and vanishes when the temperature goes beyond 
a threshold value. However, in the cases with a > a c , 
the entanglement firstly increases with temperature and 
achieves a local maximal value, then decreases with tem- 
perature and eventually vanishes. We can also find that 
in the cases with low temperature, the entanglement be- 
tween two nearest-neighbor spin ^ significantly decreases 
with a. It is natural to conjecture that the thermal state 
entanglement is closely related with fact whether this sys- 
tem has the gapped excitation. 

In Fig. 4, we plotted the log-negativity AK 1 ' 2 ) as the 
function of fc^T for three different values of a near the 
critical point. Being distinct from the entanglement of 
two spin i , if the parameter a is larger than the critical 
point, the entanglement between two nearest-neighbor 
different spin ^ and 1 decreases with the temperature 
and vanishes when the temperature goes beyond another 
threshold value which is about half of the threshold value 
concerning the two spin ^. In the case with a < a c , the 
entanglement of spin 1/2 and spin 1 firstly increases with 
the temperature, and achieves a local maximal value and 
then decreases with temperature. From Fig. 4, we can 
also find that the thermal entanglement increases with a, 
which is consistent with the ground state entanglement. 

For more deeply understanding the relation between 
thermal entanglement and the gapped excitation, the 
analysis of the entanglement of the excited state of the 
quantum mixed spin chain is necessary. In what follows, 
we briefly present a simple numerical method to calculate 
the log-negativity of the first excited state. In fact, the 
reduced density matrix of two spins in the first excited 
state can be obtained by the following procedure. 



P = 



\e-? H « ± (e-^* \g) (g\ + \ ei ){ ei \) (5) 
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FIG. 2: The log-negativity iV' 1,1 ' of the reduced density oper- 
ator of two nearest-neighbor spin i in a unit cell of the ground 
state and the log-negativity N^ 1,2 ' of two nearest-neighbor 
different SDin i and 1 are Dlotted as the function of a. The 
leng 
tern 




if the temperature is very very low. Then the reduced 
density matrix of any two spins can also be expressed 



FIG. 3: The log-negativity TV' 1 ' 1 ^ of the reduced density op- 
erator of two nearest-neighbor spin i in a unit cell of the 
thermal state is plotted as the function of fcsT. The length 
of the lattice sites is chosen as N = 128. From top to bottom, 
a = 0.758, 0.768, 0.778, respectively. 



as the weight sum of the reduced density matrices of the 
ground state and the first excited state. By iterating this 
procedure, we can obtain the log-negativity of two spins 
in any excited state. The details of this work will be 
presented elsewhere. 



V. CONCLUSSION 

In this paper, we we investigate the ground state entan- 
glement of the mixed spin chain up to 128 sites by explor- 
ing the log-negativity as an entanglement measure. Two 
distinct behaviors of pairwise entanglement are found. 
We provide a evidence that the ground state entangle- 
ment between two 1/2 spins or between 1/2 spin and 1 
spin are closely related to the so-called VBS phase tran- 
sition. At the point of VBS phase transition, the ground 
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FIG. 4: The log-negativity JV' 1 ' 2 - 1 of the reduced density oper- 
ator of two nearest-neighbor spin i and 1 in a unit cell of the 
thermal state is plotted as the function of ksT. The length 
of the lattice sites is chosen as N — 128. From top to bottom, 
a = 0.778, 0.768, 0.758, respectively. 



state between two 1/2 spins can be transferred to two 
different kinds of spins i.e. the spin 1/2 and spin 1 or 
vice versa. Furthermore, the thermal state pairwise en- 
tanglement is also studied and whether the thermal en- 
tanglement can be improved by increasing the tempera- 
ture or not is also addressed, which strongly depends on 
the parameter a = J2/ J\- It is shown that A^ 1,1 ' de- 
creases from 1 to zero as the parameter a increases from 
to 1. At the point of VBS phase transition labelled 
by a c ~ 0.768, the entanglement experience a sharply 
descent. While the entanglement between two nearest- 
neighbor different spins i and 1 characterized by log- 
negativity is zero for small values of a, and near 



the critical point it increases with a and experience a 
sharply ascent at the critical point. Unfortunately, we 
find that the ground state entanglement only exists be- 
tween two spins in the same unit cell. In the numerical 
calculation of the log-negativity of the ground state, we 
have set j3 = 200 and N — 128. if the parameter a 
is smaller than the critical point a Cl the entanglement 
decreases with the temperature and vanishes when the 
temperature goes beyond a threshold value. However, in 
the cases with a > a c , the entanglement firstly increases 
with temperature and achieves a local maximal value, 
then decreases with temperature and eventually vanishes. 
We can also find that in the cases with low temperature, 
the entanglement between two nearest-neighbor spin ^ 
significantly decreases with a. Being distinct from the 
entanglement of two spin i, if the parameter a is larger 
than the critical point, the entanglement between two 
nearest-neighbor different spin ^ and 1 decreases with 
the temperature and vanishes when the temperature goes 
beyond another threshold value which is about half of the 
threshold value concerning the two spin i. In the case 
with a < a c , the entanglement of spin 1/2 and spin 1 
firstly increases with the temperature, and achieves a lo- 
cal maximal value and then decreases with temperature. 
It is natural to conjecture that the thermal state entan- 
glement is closely related with fact whether this system 
has the gapped excitation. 
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